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Abstract 

This note gives conditions that must be imposed to algebraic multilevel discretizations involving at the same time 
nodal and edge elements so that a gradient-prolongation commutativity condition will be satisfied; this condition 
is very important, since it characterizes the gradients of coarse nodal functions in the coarse edge function space. 
They will be expressed using graph theory and they provide techniques to compute approximation bases at each 
level. To cite this article: A. Namel, A. Name2, C. R. Acad. Sci. Paris, Ser. I 340 (2005). 

Resume 

Commutativite entre gradient et prolongement et theorie des graphes Cette note donne des conditions 
qui doivent etre imposees aux discretisations multiniveau algebriques en elements finis nodaux et d'arete de facon 
a assurer la commutativite entre gradient et prolongement ; cette relation importante caracterise les gradients 
des fonctions nodales grossieres dans l'espace des fonctions d'arete grossieres. Ces conditions seront exprimees en 
terme de graphes et elles permettent d'introduire des methodes de calcul des bases d'approximation aux differents 
niveaux. Pour citer cet article : A. Namel, A. Name2, C. R. Acad. Sci. Paris, Ser. I 340 (2005). 



Version francaise abregee 

L'approximation numerique du champ electrique ou magnetique utilise frequemment les elements finis 
d'arete dont la relation avec les elements finis nodaux traduit des proprietes importantes au niveau 
discret [1]. Dans ce qui suit, nous considererons les elements de plus bas degre : Pi en nodal et ordre 1 
incomplet pour les aretes. Des qu'on traite des problemes de grande taille, une strategie multiniveau est 
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un choix interessant. Pour les systemes provenant de discretisations par elements finis d'arete, Hiptmair 
a introduit des methodes multiniveau pour une hierarchie de maillages emboites [2]. 

Cependant, dans des applications realistes, on ne dispose generalement pas de maillages structures. La 
strategic multiniveau algebrique va done s'imposer : il s'agit de definir des fonctions grossieres nodales et 
d'arete grace aux contributions de paquets de fonctions fines nodales et d'arete ; les combinaisons lineaires 
(la) et (lb) defmissent respectivement ces fonctions grossieres nodales et d'arete. 

Par construction les gradients des fonctions nodales fines appartiennent a l'espace des fonctions d'arete 
fines ce que traduit la relation (2). Dans cette relation, G h est la matrice d'incidence arcs-sommcts du 
graphc oricntc naturellement associe au maillage de travail. Les orientations des arcs sont arbitraires. 

Pour adapter aux methodes algebriques les lisseurs des methodes geometriques, Reitzinger et Schoberl [3] 
ont introduit une representation explicite des gradients des fonctions grossieres nodales dans la base des 
fonctions grossieres d'arete, donnee par la relation (3) ou G H est une matrice d'incidence arcs-sommets. 

En regroupant les relations (1) a (3), nous obtenons la relation matriciellc (4). La matrice a est 
construitc par cxemplc par les methodes definies dans [4] qui pcrmettent d'obtenir les fonctions grossieres 
nodales commc partition de l'unite et de contraindre leurs supports a etre inclus dans des ensembles 
geometriques convenablement choisis. 

Connaissant G h et a, nous souhaitons choisir G H comme matrice d'incidence arcs-sommets d'un graphc 
oriente S H . Nous donnons dans cette note une condition necessaire et suffisante sur ce graphe, la propo- 
sition (2.3), qui assure l'existence d'une solution de (4). En effet, nous associerons, par un precede decrit 
dans la partie en anglais, a chaque arete fine un sous-graphe du graphc grossier, qui doit etre connexe. 

La connaissance de ces sous-graphes donne les degres de libcrte disponibles pour determiner des fonc- 
tions d'arete grossieres compatibles avec les fonctions nodales grossieres ; en resolvant un probleme de flot 
sur ces sous-graphes, voir (14), nous pouvons alors construire la matrice (3 (Section 4). 



1. Introduction 

Numerical approximation of electric or magnetic field uses often edge finite elements whose relation with 
nodal finite elements contains important properties at discrete level [1]. In this note we restrict ourselves 
to lowest order approximation : P\ for nodal elements and incomplete order 1 for edge elements. In order 
to solve large problems, multilevel methods are an attractive choice. While, for systems coming from 
edge element discretisation, Hiptmair [2] proposed multilevel methods using nested meshes, engineering 
applications do not usually provide structured meshes. Therefore, algebraic multilevel methods are an 
interesting option: we have to build coarse nodal and edge functions by using aggregates of fine nodal 
and edge functions. If {(j)p) p= i M h an d {^i)i=i,....E h respectively denote fine nodal and edge bases, the 
following linear combinations define coarse nodal and edge functions: 



By construction, the gradients of fine nodal functions belong to the space of fine edge functions: 




(la) 



P =i 




(lb) 




(2) 
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where G h is the edge-node incidence matrix of the digraph naturally associated with the initial mesh. 
The orientation of the edges can be arbitrarily chosen. 

In [3] , Rcitzinger and Schobcrl deduced their smoother from the matrix G H involved in the relation: 

E H 

Vn G {1, . . . , N H }, grad(^) = £ Gf„A? , (3) 

e=l 

which states that the gradients of the coarse nodal functions must also belong to the space of coarse edge 
functions. The matrix G H is an edge-node incidence matrix as in the structured case. Relation (3) does 
not guarantee the efficacy of the algebraic multilevel method but it leads to relevant strategies. 
Gathering Equations (1), (2) and (3), we obtain the matrix relation: 

G h a = [3G H . (4) 

The matrix a is constructed following for instance the methods defined in [4] , which provides a family of 
coarse nodal functions, making up a partition of unity, whose supports satisfy appropriate conditions. 

Knowing the left-hand side of (4) , we want to choose G H as an edge- node incidence matrix of a digraph 
S H , and we will give conditions on the coarse graph S H , which ensure the existence of a matrix satisfying 
(4). Moreover, the proof of the proposition indicates how to choose the degrees of freedom which enables 
us to define the coarse edge functions. It also helps us to construct 0. 



2. Notation and statement of the problem 

Let {L n ) n =\,...,N H be sets of indices in {1, ... , N h } such that: 

N H 

\jL n = {l,...,N h }. (5) 



The matrix a describes the coarse nodal basis; we assume that it is has been previously computed and 
it has the following properties: 

- the coarse nodal functions make up a partition of unity, which can be algebraically stated as: 



N 



H 



Vpe{l,...,N h }, ^a pn = l, (6) 

n=l 

- in order to restrict the support of each coarse basis function <j>^ , the indices of the non-zero components 
of (f)H are included in the set L n , i.e.: 

P e{l,...,N h }\L n => a pn =0. (7) 

The fine nodal function <^ contributes to the coarse nodal function (f>^ if p belongs to L n . 

We have a reciprocal set- valued function L: the set L p is the set of coarse nodal function indices to which 
the fine nodal function (j^ contributes. For the fine graph in Figure 1(a), we set L\ = {1,2,3,4,5,6,7}, 

L2 = {5, 6, 8, 9, 13, 14} and L3 = {7, 8, 10, 11, 12}. One obtains, for instance, the set L7 = {1, 3}. 

We define two families of sets of fine edge function indices. We will denote a directed fine edge i by pq h 
where p and q are respectively the starting and ending nodes of the edge i. A similar notation is used for 
a directed coarse edge e = mn H . 

The set C n is the set of indices of fine edges which have an extremity in L n : 

C n = {i£{l,...,E h }:i=pq h , peL n orqeL n }. (8) 



The fine edge function X^ contributes to the gradient of the coarse nodal function if i belongs to 
C n . Indeed, for the directed fine edge i — pq h , G\ r is equal to —1 if r =p and +1 if r = q. Moreover, if p 
and q are not in L„, the components a pn and a qn vanish according to (7); therefore: 

i G {1, . . • , E h } \ C n (G h a. n )i = 0, (9) 

where a, n denotes the n-th column of a. The reciprocal set-valued function C is such that Ci is the set 
of coarse nodal function indices to whose gradient the fine edge function X^ contributes. On Figure 1(b), 
the fine edges are numbered, set C3 is highlighted and we can note, for instance, the set C$ = {1,3}. 
Let e = rrvn H be an edge of the coarse graph S H ; we define: 

h = c n n c m . (10) 

By analogy with the structured case and for restricting the support of Ag , we enforce: 

ie{l,...,E h }\l e =► ft e = 0. (11) 

The fine edge function X^ contributes to the coarse edge function X^f if i belongs to I e . The set-valued 
function I is such that J, is the set of coarse edge function indices to which the fine edge function X^ 
contributes. The coarse graph in Figure 1(c) is related to the fine in Figure 1(a). Set I e3 is represented in 
Figure 1(d). 




rounded. edges. 



Figure 1. Representation of the fine and coarse graphes, sets (in)n=l,...,3i C3 and I e3 . 

The following statement can be easily deduced from (8) and the definition of G h : 
Lemma 2.1 If i denotes the edge pq , Ci — L p U L q . 

In order to simplify notations, we introduce the set F = \i G {1, . . . , E h } : 1^ ^ 0}, since some fine 
edge functions might not contribute to any coarse edge functions. 

For any fine edge i, let S H ' 1 be the induced subgraph defined by Q: the vertices of S H l are the vertices 
of S H , which are indexed by the elements of Ci and the edges of S 11,1 are those edges of S H whose 
extremities are vertices of S H ' 1 . 

The following lemma is a direct consequence of definition (10): 
Lemma 2.2 For any edge i G F, the edges of S '* are those edges of S H which are indexed by I{. 

We may now state precisely our main result, which gives a necessary and sufficient condition on the 
coarse graph S H permitting the resolution of (4): 

Proposition 2.3 For all matrices a satisfying conditions (6) and (7), there exists a matrix [3 satisfying 
(11) and solving (4) iff for all i, the induced subgraph S H l is connected. 
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3. The essential steps of the proof 



First step. Many relations in (4) reduce to = 0: this is the case for n £ Ci. 

Indeed, according to (9) the (i,n) coefficient of the right-hand side of (4) vanishes. 

Conversely, if e does not belong to 7j, according to (11) and the definition of Ii, /?j e vanishes and: 

E H 

£AeGf n = ^AeGf„. (12) 

On the other handjf the directed coarse edge e denoted by Im belongs to 7j, Lemma 2.2 implies that I 
and m belongs to Ci. However, for G^ n not to vanish for all e, m or I must be equal to n, which means 
that n belong to d, and this contradicts the assumption n ^ d. 

Second step. We look at all the other equations, i.e. those for which n e Cj. We note that (12) remains 
and that the edges indexed by Ii are precisely those of the graph S H ' 1 according to Lemma 2.2. 

We assume now i £ f and we define G H l as the edge-node incidence matrix of S H l and the (i,n) 
equation of (4) is rewritten: 

E AeGf/ = Qi, n where 9 4 ,„ = E G h „a rn . (13) 

eeh reL " 
This could be satisfied for all couples (i, n) such that n 6 {1, . . . , N H } and i e C n or equivalcntly 
i E {1, . . . , E h } and n e Ci. For a fixed i, we may write that /3, # , the i-th row of (3 satisfies the system: 

E AeGf„ l = ©,,„, Vn G Ci. (14) 
eeZ 

Thus, we solve line by line for (3 and we see that (14) is a flow problem whose solution is of the form: 

(15) 

with (f3" m Y G ker(G ff ' 1 )' and a particular solution. 

More precisely, let T* be a spanning tree for S H ' Z ; call the edge- node incidence matrix associated 
with T l ; we know that T 4 has \Ci\ — l rows and \Ci\ columns, and it is of rank \Ci\ — 1. We choose a vertex 
m in T l and we solve the system: 

E 4^ = e iin , Vn G a \ {m}, (16) 

where £(T l ) denotes the set of indices of the edges of T % . The system (16) is a regular system of \C%\ — 1 
equations with |Cj| — 1 unknowns, and we put (3' ie equal to if e is in Ii \ £(T l ). 

It remains to show that the forgotten equation of index m in (16) is automatically satisfied. Indeed, by 
denoting i by pq h , we sum the right-hand side of (14) with respect to n G Cf. 

E ©J,™ = e a 9« - = °' ( 17 ) 

neCi n£L p UL q 

since in view of (6) and (7), E a pn = E a ? n = ^ ano - E a <?™ = E a< ?" = ^' 

n£L p UL q n£L p n£L p UL q n£L q 

On the other hand, if we sum the left-hand side of (14) with respect to n G Ci, we obtain: 

£ £ = £ £ G - = 0, (18) 
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since each line of G H ' 1 contains only two non-zero coefficients +1 and — 1. 

E H 

If i £ F, \Ci\= \L p \= \L q \= 1 and the relation \_. PieG^ n = &i >n is satisfied from (12) and (17). 

e=l 

Now we assume that S 11 ' 1 is not connected and we denote by the nodes of a connected component. 

E H 

For the same reasons as in (18), if (3 satisfies (11) one gets f3 ie G^ n = 0. 



n6C, 6=1 



However we can construct a matrix a satisfying (6) and (7) such that G^ r a rn ^ 0. In fact, in 

view of (7) , for i = pq h we can write: 

^ ^ ^ ^ G \ r (y. rn ^ ^ Qqn Qpn ^ ^ ^qn ^ ^ ^pn ■ (1^) 



neCi' -1 nee, c,nL„ CinL„ 



Since Ci is strictly included in L p U L q , we will have Ci C\ L p ^ L p or Ci C\ L q ^ L 9 . Depending on the 
situation, we can construct a suitable matrix a such that: 

^ a qn = 1 and ^ a pn = 0^ or ^ ^ a 9 „ = and ^ a pn = 1 

'CiOL, CiHLp CiDLq CinL p 

For these matrices a, the condition defined by (4) cannot be ensured. 



4. Construction of the coarse edge functions 

For a coarse graph satisfying the condition of Proposition 2.3 and by using the decomposition (15), any 
compatible matrice can be written f3 = /?' + (3" , where the complete matrices are defined by gathering 
the lines of index i f} im , /?-. and (3",. The computation of each (3' i% can be done by solving system (16). 
As concerns f3" % , a basis of the kernel of (G 11 ' 1 ) 1 is given by a set of ki independent cycles of S ,% , Then, 
Y]. ~ ki degrees of freedom should be determined by minimising an appropriate energy functional; such 
a problem is introduced in [5] and can be related to explanations in [4] . 

We thank Michelle Schatzman for many fertile discussions. 
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